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Abstract. We discuss how measurements of the absorption of 7-rays from GeV to TeV 
energies via pair production on the extragalactic background light (EBL) can probe 
important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy 
formation, set within the hierarchical structure formation scenario, to obtain predic- 
tions of the EBL from 0.1 to 1000/^m. SAMs incorporate simplified physical treatments 
of the key processes of galaxy formation — including gravitational collapse and merging 
of dark matter halos, gas cooling and dissipation, star formation, supernova feedback 
and metal production — and have been shown to reproduce key observations at low 
and high redshift. Here we also introduce improved modelling of the spectral energy 
distributions in the mid-to-far-IR arising from emission by dust grains. Assuming a 
flat ACDM cosmology with fi m = 0.3 and Hubble parameter h = 0.65, we investigate 
the consequences of variations in input assumptions such as the stellar initial mass 
function (IMF) and the efficiency of converting cold gas into stars. We conclude that 
observational studies of the absorption of 7-rays with energies from ~10 Gev to ~10 
TeV will help to determine the EBL, and also help to explain its origin by constraining 
some of the most uncertain features of galaxy formation theory, including the IMF, 
the history of star formation, and the reprocessing of light by dust. 



I INTRODUCTION 

The extragalactic background light (EBL) represents all the light that has been 
emitted by galaxies over the entire history of the universe. The EBL that we 
observe today is an admixture of light from different epochs, its spectral energy 
distribution (SED) distorted by the redshifting of photons as they travel to us from 
sources at different distances. It is therefore a constraint on both the intrinsic 
SEDs of the sources and their distribution in redshift. At present, there is more 
than a factor of two uncertainty in the amplitude of the EBL in the UV, optical, 
and near- infrared [1]. The EBL in the mid-IR is even more uncertain. The far-IR 



background measured at £ lOOyum [2-5] represents at least half of the total energy 
in the EBL, yet the sources that produced it remain uncertain. 

High energy 7-ray astronomy promises to help resolve these uncertainties by 
providing independent constraints on the EBL, in the mid-IR with E 1 in the ~ 10 
TeV energy range, and in the 0.1-3 /im range with E 1 ~ 100 GeV via the new low- 
threshold instruments that will soon be available. High energy 7-rays from sources 
at cosmological distances are absorbed via electron-positron pair production on the 
diffuse background of photons that comprises the EBL. Thus, 7-ray observations 
of objects with known redshift and intrinsic spectral shape will constrain the EBL 
in these crucial wavelength regimes by measuring the optical depth of the Universe 
to photons of various energies. This in turn will help to constrain some of the most 
fundamental uncertainties in physical models of galaxy formation. 

In order to illustrate this, in this paper we use a "forward evolution" approach, 
which attempts to model the essential features of galaxy formation using simple 
recipes. These semi-analytic models are set within the modern Cold Dark Matter 
(CDM) paradigm of hierarchical structure formation, and trace the gravitational 
collapse and merging of dark matter halos, the cooling and shock heating of gas, 
star formation, supernovae feedback, metal production, the evolution of stellar pop- 
ulations and the absorption and re-emission of starlight by dust. This machinery 
has been used extensively to predict optical properties of low-redshift galaxies, with 
good results (e.g., [6,7]; reviewed and extended in [8,9], hereafter SP and SPF). A 
semi-analytic approach was also used by Devriendt and Guiderdoni [10] to make 
predictions of counts and backgrounds in the mid-to-far-IR, with more detailed 
modelling of dust extinction and emission, but less detailed modelling of merging 
and star formation. We have now combined the strengths of these two models, by 
integrating the stellar SEDs and dust modelling of [11,10] into the galaxy formation 
code of the Santa Cruz group. 

Some parts of the "standard paradigm" of galaxy formation represented by our 
SAMs are relatively solid. For example, once a cosmological model and power 
spectrum are specified, it is straightforward to compute the gravitational collapse 
of dark matter into bound halos using A^-body techniques, and analytic formalisms 
such as those used in our modelling [12] have been checked against these results 
[13]. Within the range of values for the cosmological parameters allowed by existing 
observational constraints (i.e., fi mattcr ~ 0.3 — 0.5, fi mattcr + Q\ ~ 1, H ~ 60 — 80 
km/s/Mpc; see e.g. [14] for a summary), these results do not change significantly. 
Similarly, modelling of gas cooling appears to be fairly robust and agrees well with 
hydrodynamic simulations [15]. However, other aspects, notably the efficiency of 
conversion of cold gas into stars, the effect of subsequent feedback due to supernovae 
winds or ionizing photons, the stellar initial mass function (IMF), and the effects 
of dust, remain highly uncertain, and some predictions are quite sensitive to their 
details. 

For example, SPF showed that the star formation history of the Universe and 
the number density of high redshift z ^ 2 "Lyman-break" galaxies (LBGs; e.g. [16]) 
may be quite different depending on whether star formation is primarily regulated 



by internal properties, such as gas surface density in a quiescent disk, or triggered 
by an external event such as an interaction. Because the largest samples of LBGs 
are primarily identified in the rest UV, model predictions are also quite sensitive 
to the high-stellar-mass slope of the IMF, and to dust extinction. At the other 
end of the spectrum is the sub-mm population detected by SCUBA, believed to be 
predominantly high redshift (z £ 2) luminous and ultraluminous infrared galaxies 
(LIRGs and ULIRGs) powered by star formation rates of hundreds to thousands 
of solar masses per year (e.g., [17]). Theoretical predictions of the numbers and 
nature of these objects are highly sensitive to the same issues (the dominant mode 
of star formation, dust, the IMF), but provide a crucial counter-balance to the 
optical observations. However, the current mismatch between the sensitivity and 
spatial resolution of optical and sub-mm instrumentation has made it difficult to 
establish the connection between the two populations observationally. 

The Milky Way, like most nearby galaxies, emits the majority of its light in optical 
and near-IR wavelengths; only about 30% of the bolometric luminosity locally is 
released in the far-infrared [18]. This was generally believed to be typical of most 
of the starlight at all redshifts until the discovery of the far-IR part of the EBL 
by the DIRBE and FIRAS instruments on the COBE satellite, at a level ten times 
higher than the no-evolution predictions based on the local luminosity function of 
IRAS galaxies, and representing twice as much energy as the optical background 
obtained from counts of resolved galaxies [19]. This result suggests that either the 
dust extinction properties of "normal" galaxies change dramatically with redshift, 
or a population of heavily extinguished galaxies (perhaps analogous to local LIRGs 
and ULIRGs) is much more common at high redshift than locally, or both. Some 
of these galaxies may have already been observed, at 15 /im by ISO [20], and at 
850 /im by SCUBA [21]. 

Guiderdoni et al. [22,10] showed that their simplified semi-analytic model could 
reproduce the multi-wavelength data only if they introduced a population of heav- 
ily extinguished galaxies with high star formation rates, and with strong evolution 
of number density with redshift. This population was introduced ad-hoc by [22,10], 
but as discussed by these authors, by [23] (based on [24]), and also by SPF, the 
increasing importance of starbursts at high redshift, due to the increasing merger 
rate and higher gas fractions, is a natural mechanism to produce this population. 
The models of SPF contain a detailed treatment of mergers and the ensuing colli- 
sional starbursts, which has been calibrated against the merger rate in cosmological 
./V-body simulations [25] and the starburst efficiency in hydrodynamical simulations 
[26,27]. Moreover, they produced good agreement with observations of LBGs and 
Damped Lyman-a systems (SPF) as well as low redshift galaxies (SP). Therefore, 
it will be extremely interesting to see if these same models, when combined with 
the more sophisticated treatment of dust extinction and emission developed by De- 
vriendt, Guiderdoni, and collaborators, will be able to simultaneously reproduce 
observations over the broad range of wavelengths and redshifts discussed above. 

In the next section we briefly describe the ingredients of our models, and then 
present the results of the predicted EBL. The following section presents the impli- 




FIGURE 1. (a) The star formation rate density predicted by our models, for two different 
recipes of star formation. Both models produce about the same total mass density of stars by 
z — (i.e., the areas under the curves are equal when they are plotted linearly vs. time), but the 
collisional starburst model (CSB) peaks at higher redshift. (b) Comoving luminosity density at 
2000 A as a function of redshift. Data points represent the observed global luminosity density at 
rest ~ 2000A, obtained by integrating the observational best-fit Schechter luminosity functions 
over all luminosities (pl = 0*£*r(2 — a)). The z — 0.15 point is from [28], the z ~ 0.4 and 1.2 
points are from [29], and the z ~ 3 and z ~ 4 points are from [16]. The curves for our four models 
arc labeled as in Figure 3. The model curves have been corrected for dust extinction using the 
approach described in the text. 

cations for 7-ray attenuation, and the final one briefly discusses some alternative 
treatments and our own conclusions. The work summarized here is a brief, prelim- 
inary sample of the results which will soon be presented in a series of papers, now 
in preparation, on the EBL and its breakdown into various kinds of sources and on 
the implications for 7-ray astronomy. 

II SEMI-ANALYTIC MODELLING 

In this section we briefly describe the ingredients of our models. Readers can 
refer to SP and SPF for more details, and to [30] for a brief introduction. 

Using the method described in [12], we create Monte-Carlo realizations of the 
masses of progenitor halos and the redshifts at which they merge to form a larger 
halo. These "merger trees" (each branch in the tree represents a halo merging 
event) reflect the collapse and merging of dark matter halos within a specific cos- 
mology We truncate the trees at halos with a minimum circular velocity of 40 
km/s, below which we assume that the gas is prevented from collapsing and cool- 
ing by photoionization. Each halo at the top level of the hierarchy is assumed to 
be filled with hot gas, which cools radiatively and collapses to form a gaseous disk. 



The cooling rate is calculated from the density, metallicity, and temperature of the 
gas. Cold gas is turned into stars using a simple recipe, depending on the mass of 
cold gas present and the dynamical time of the disk. Supernovae inject energy into 
the cold gas and may expell it from the disk and/or halo if this energy is larger 
than the escape velocity of the system. Chemical evolution is traced assuming a 
constant yield of metals per unit mass of new stars formed. Metals are initially 
deposited into the cold gas, and may later be redistributed by supernovae feedback, 
and mixed with the hot gas or the diffuse (extra-halo) IGM. 

When halos merge, the galaxies contained in each progenitor halo retain their 
seperate identities until they either spiral to the center of the halo due to dynamical 
friction and merge with the central galaxy, or until they experience a binding merger 
with another satellite galaxy orbiting within the same halo. All newly cooled gas is 
assumed to initally collapse to form a disk, and major (nearly equal mass) mergers 
result in the formation of a spheroid. New gas accretion and star formation may 
later form a new disk, resulting in a variety of bulge-to-disk ratios at late times. 

For an assumed IMF, the stellar SED of each galaxy is then obtained using stellar 
population models. Here we use the multi-metallicity stellar SEDs of [11] for the 
Salpeter and Kennicutt IMF cases, and the solar metallicity GISSEL models [34] 
for the Scalo IMF. (We have found that using evolving metallicity rather than 
solar metallicity SEDs has a relatively small impact on the resulting EBL.) Dust 
extinction is modelled using an approach similar to that of [10]. The optical depth 
of the disk is assumed to be proportional to the column density of metals. We 
then use a simple slab geometry where stars and gas are homogenously mixed, and 
assign a random inclination to each galaxy to compute the absorption. We use a 
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FIGURE 2. The three stellar Initial Mass Functions (IMFs) used here: Kennicutt [40], Salpeter 
[39], and Scalo [38]. 
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FIGURE 3. Comoving luminosity density at redshift z = as a function of wavelength. Data 
points represent the observed global luminosity density of the Universe, obtained by integrating 
the observational best-fit Schcchtcr luminosity functions over all luminosities (pl — </>*£*r(2 — a)). 
The far-UV point (at A = 200(M) is from the luminosity function from FOCA observations [28], 
extrapolated from the mean redshift of the sample (z = 0.15) to z = assuming that the 
luminosity density scales with redshift as pl oc (1 + z) , as indicated by the observations of 
Cowie et al. [29]. The B-band point (A = 4400A) is from the luminosity function derived from 
the 2dF redshift survey [31], the R-band point (A = 6940A) is from the Century Redshift Survey 
[32], and the K-band point (A = 2.2pm) is from [33]. The model curves are obtained by simply 
summing the spectra of all z = galaxies in our models with the appropriate Press-Schechter 
weighting. 

metallicity dependent extinction curve, following [22,10]. 

All absorbed light is re-radiated at longer wavelength. In general, any galactic 
dust emission spectrum can by represented by a combination of three components: 
1) hot dust (as in Hn regions), 2) warm dust (as in the diffuse Hi), and 3) cold dust 
(as in molecular clouds). In the models of Devriendt et al. [11], these components 
are modelled as a mixture of polycyclic aromatic hydrocarbon molecules (PAH), 
very small grains, and big grains. Big grains may be either cold (~ 17 K), or heated 
by radiation from star-forming regions (as suggested by observations of typical 
local starburst galaxies like M82). A set of template spectra is then constructed 
for galaxies of varying IR luminosity, with admixtures of the various components 
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selected in order to reproduce the observed relations between IR/sub-mm color and 
IR luminosity. A similar approach was used by [35], using a mixture of a typical 
Orion-like Hu spectrum and an Hi spectrum constructed to fit DIRBE observations 
of the diffuse ISM [36]. Here, we use the more empirical emission templates of [35] 
(kindly provided in electronic form by E. Dwek), but we obtain very similar results 
with the models of [11]. 

The recipes for star formation, feedback, chemical evolution, and dust optical 
depth contain free parameters, which we set for each model (see SP) by requiring 
an average fiducial "Milky Way" galaxy to have a K-band magnitude, cold gas 
mass, metallicity, and average B-band extinction as dictated by observations of 
nearby galaxies. 

Figure la shows the global star formation rate density for the two star formation 
recipes that we consider here. The "fiducial" model is the collisional starburst 
(CSB) model favored by SPF, in which bursts of star formation may be triggered 
by galaxy collisions. The "Late star formation" model is the Constant Efficiency 
Quiescent (CEQ) model of SPF, in which cold gas is converted to stars only in a 
quiescent mode with constant efficiency This produces a star formation history 
similar to the models of [37], in which the peak in the star formation history occurs 
at a considerably later epoch (z ~ 1.5) than in the CSB model. Figure 3 shows 
the resulting luminosity density as a function of wavelength at z — 0. For the CSB 
model, we consider three different choices of IMF: Scalo [38], Salpeter [39], and 
Kennicutt [40]. These IMFs are graphed in Figure 2. For the CEQ model we show 
only the Kennicutt case. This is compared with the observed luminosity density 
from nearby galaxies, obtained by integrating the luminosity functions of galaxies 
resolved in recent redshift surveys at wavelengths ranging from 0.2 to 2.2 jum. The 
spikes in the model predictions at ~ 5 — 12/xm are caused by the PAH features 
mentioned above. All four models, when normalized to the observed K-band Tully- 
Fisher relation, produce reasonable agreement with the observed luminosity density 
in the B and K bands. 1 This is perhaps not surprising, yet it was not guaranteed. 
However, there is a noticable difference in the far-UV and the mid- to far-IR. 
The Scalo IMF produces too little UV light relative to optical and near-IR light, 
whereas the Kennicutt and Salpeter IMFs are in much better agreement with the 
data. These IMFs produce more high mass stars than the Scalo IMF, and thus 
more ultraviolet light to be absorbed and re-radiated by dust in the far IR. In 
Fig. lb we show the redshift evolution of the far-UV (2000A) luminosity density, 
compared with observations. The Scalo model falls short at all redshifts, and the 
CEQ model, which agrees at z — 0, falls short at higher redshifts. It is encouraging 
that our very simple model for dust extinction, which we normalized in the B-band 
at z — 0, appears to yield the appropriate level of dust extinction in the UV at 

^ In [30], we renormalized all the models by requiring that they all agreed with the K-band point 
at 2.2 /urn. Here we do not do this since, as Fig. 3 shows, the SAM parameters chosen for each 
case to produce an average fiducial "Milky Way" as described above are already in agreement 
with this data within the errors. Also, our current SAMs [8,9] use a corrected version [41] of the 
Press-Schechter formalism, which obviates our previous motivation for the K-band normalization. 



higher redshifts (SPF). 



Ill THE INTEGRATED EXTRAGALACTIC 
BACKGROUND LIGHT 

Figure 4 shows the EBL produced by our four models, obtained by integrating 
the light over redshift (out to z — 4) with the appropriate K-corrections due to 
cosmological redshifting. We compare this with a compilation of observational 
limits and measurements of the EBL. It is apparent that there is at least as much 
energy in the far-IR part of the EBL as in the entire optical and near-IR bands. For 
example, Puget and collaborators [1] estimated that the total energy in the EBL is 
between 60 and 93 nW m~ 2 sr^ 1 , with between 20 and 41 nW m~ 2 sr _1 contributed 
by the optical and near-IR, and between 40 and 52 nW m~ 2 sr _1 coming from the 
far-IR. If the possible detection of the EBL at 60 [xm by Finkbeiner et al. [48] is 
correct, that would further increase the far-IR EBL; however, as Puget discussed 
in his talk at this conference, it is very difficult to determine the EBL at 60 /xm 
since the zodiacal light is so much brighter at that wavelength. 

In units of critical density p c , Qebl = (47r/c) (Jebl / PcC 2 ) = 2.5 x 10~ 8 /ebl^~ 2 , 
where Jebl is in units of nW m~ 2 sr _1 . The total energy density in the EBL 
corresponding to the lower and upper estimates of [1] is ^ebl = (3.6 — 5.5) x 
10~ 6 (/i/0.65)~ 2 . Although the EBL includes energy radiated by active galactic 
nuclei (AGNs) as well as stars, it is unlikely that AGNs contributed more than a 
few percent of the total. The total energy radiated by AGNs is = VPbrc 2 , 

where the efficiency of conversion of mass to radiated energy in AGNs is rj ~ 
0.05. Correspondingly, = ^bh(1 + ^bh) _1 ~ 4.5 x 10- 8 h-\r]/0.05)[3/(l + 

zbe)] ~ 0.02f2EBL- 2 So for simplicity, in this paper we will neglect the contribution 
of AGNs to the EBL. 

Several interesting features emerge from the comparison of our SAM models 
with the EBL data. In the UV to near-IR, the models are much closer to the direct 
measures of the EBL obtained by [44-46], although the Scalo IMF produces less 
light in the UV because it has fewer high-mass stars. Recall that the Kennicutt 
model agreed well with the observed luminosity density at z — 0, and the observed 
redshift evolution of the luminosity density in the rest UV. This suggests that the 
extra factor of 2-3 in the direct measurements of the EBL must arise from a rapidly 
evolving population of objects which are too faint or too low in surface brightness 
to be detected in the samples used to obtain the counts (e.g., [19]). We are in 
the process of attempting to determine whether the observational selection effects 
inherent in the measured counts are sufficient to explain this discrepancy for our 
modelled population. A second interesting point is that all the models satisfy the 
lower limits from counts in the mid-IR (15 /iin; [20]) and the sub-mm (850 /im; 

2) Updating [49], we have estimated fi B H = (MBH/M S p hcro i d )fi S p hcroid « (1.5 x 1(T 3 )(1.8 x 
10 _3 /i _1 ), using the observed (loose) correlation [50] between a black hole mass and that of the 
galactic spheroid in which it is found, and the estimated cosmological density of spheroids [51]. 
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FIGURE 4. Extragalactic background light: models and data. The far-UV points are from 
STIS (inverted filled triangles) [42] and FOCA observations (filled triangle) [43]. The lower optical 
points (filled squares) arc lower limits from resolved sources [19]; the upper ones (open diamonds) 
are from absolute photometry [44]. The near-IR points are from DIRBE: (open circle) [45], (open 
triangles) [46]. The point at 15 /mi is from ISOCAM resolved sources [20], and is thus a lower 
limit. The far-IR points are from DIRBE (filled circles) [4,47], (stars) [48]. The curves are our 
results from modelling the history of star formation in the ACDM cosmology using semi-analytic 
methods: a model with both quiescent star formation with constant efficienty and starbursts, 
with Kennicutt, Salpeter, and Scalo IMFs, and a model with only quiescent star formation with 
constant efficienty (Late SF). The lower light dotted curve is the ACDM EBL calculated using 
our previous methods [30] for the Salpeter IMF, and the upper one is the same curve to 80 /im 
multiplied by 2.5 for comparison with Mrk 501 data as analyzed by [53] (see text). Note that 
10~ 6 erg s _1 cm~ 2 sr _1 = 1 nW m~ 2 sr _1 . 



[21]). Of our four new EBL curves, the Late SF model and the fiducial Kennicutt 
model are also consistent with the DIRBE/FIRAS measurements at 140 and 240 
/im. The models differ significantly in the mid-IR, ~ 10 — 60yum, where the EBL 
can be probed by TeV 7-rays. The lower dotted curve, representing our previous 
attempt [30] to model the EBL, is well below the 15 /xm lower limit as well as the 
DIRBE measurements at longer wavelengths. As we stated in [30] , we expected our 
EBL results to change as we improved our dust emission modelling. In addition to 
inclusion of the PAH features, the new dust emission model has more warm dust 
than the one used in [30]. 

We now discuss constraints from the TeV 7-ray observations. 



IV ATTENUATION OF HIGH-ENERGY 7-RAYS 

Figure 5 shows the 7-ray attenuation predicted by the four ACDM models con- 
sidered here, for sources at redshifts z s = 0.03 and 0.10. All of the models predict 
rather little absorption at £ 5 TeV for sources at z s = 0.03, but fairly sharp cut- 
offs above 10 TeV, especially for the Late SF model. That model may be in conflict 
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FIGURE 5. The attenuation factor, exp(— r) for 7-rays as a function of 7-ray energy for the 
four ACDM models considered in Fig. 4. The assumed redshift of the source, z s , is indicated for 
each set of curves. 



with the data from Mrk 501. The synchrotron self-Compton (SSC) model, in which 
~keV synchrotron X-radiation from a very energetic electron beam is Compton up- 
scattered by the same electrons to produce the observed ~TeV 7-rays, appears to 
explain both the keV-TeV spectra and their time variation for the blazars Mrk 421 
and 501, both at z ~ 0.03 (see, e.g., [53,54] and references therein). Using a sim- 
plified SSC model and keV X-ray data to predict the unattenuated TeV spectrum 
of Mrk 501, Guy et al. [53] used CAT and HEGRA data to estimate the amount 
of 7-ray attenuation. They find that there is a rather good fit to the observed 
attenuation for the ACDM-Salpeter EBL from [30] when it is scaled upward by a 
factor of up to about 2.5 across the wavelength range 1-80 jum; this is the upper 
Salpeter curve on Fig. 4. The [53] 1-er upper limit for 20-80 /xm is a scaling factor 
of 3.4. Our new Salpeter curve appears to be rather consistent with this rescaling 
of our old Salpeter one, the Kennicutt curve may be a little high, but the Late SF 
curve appears to be definitely too high. 

The compatibility of our new EBL calculations with the available data on TeV 
7-ray attenuation is definitely worth further investigation. The results appear to 
be sensitive to the details of the models, raising the hope that they may be able to 
help answer important questions about star formation and dust reradiation, and 
also help to test the SSC modelling. For example, Figure 6 shows the optical depth 
as a function of 7-ray energy E 1 , compared with 7-ray attenuation results from Mrk 
501 with detailed SSC models. This figure is like Fig. 10 of Krawczynski et al. [54]. 
While that figure showed OD(£' T ) - OD(0.5 TeV), following Krawczynski's advice 
we here plot OD(£' T ) - OD(1.5 TeV), normalizing to the data at 1.5 TeV since the 
systematic error in the curvature of the spectrum strongly increases below 1.5 TeV, 
corresponding at 0.5 TeV to a flux uncertainty of 50%. (Krawczynski also kindly 
updated his model curves for this figure to take into account the 15% HEGRA 
energy uncertainty. We followed his advice to omit the highest-energy point, which 
has a statistical significance well below 2a.) The conclusions from Fig. 6 appear 
to be consistent with those from our discussion of [53]: the higher model curve 
seems compatible with our new Salpeter results, taking into account that the error 
bars on the data points also apply to the model curves; the Kennicutt model also 
appears to be consistent, except perhaps at the highest E 7 ; and the Late SF model 
definitely predicts too much attenuation. 

Figure 7 depicts the 7-ray "absorption edge," the redshift of a source corre- 
sponding to an optical depth of unity, as a function of 7-ray energy. Travelling 
through the evolving extragalactic radiation field, 7-rays from sources at lower red- 
shift suffer little attenution. The universe becomes increasingly transparent as E 1 
decreases, probing the background light at increasingly short wavelengths. (We are 
using the treatment of [52] to account for absorption of ionizing radiation by the 
Lyman alpha forest.) The models all have the same qualitative features, but differ 
significantly quantitatively. The location of the absorption edge is affected both by 
the assumed IMF and by the history of star formation. There is more absorption 
at most redshifts with the Kennicutt IMF because with a higher fraction of high 
mass stars, it is more efficient at producing radiation for a given stellar mass; there 
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FIGURE 6. Increase of the optical depth (OD) for 7-rays due to intergalactic extinction, inferred 
from comparison of the observed Mrk 501 spectrum from HEGRA with that estimated using SSC 
models. The points correspond to Doppler factor Sj = 25, magnetic field B = 0.037 G; the upper 
and lower light solid curves correspond respectively to (Sj,B) — (100,0.012 G) and (25,0.12 G), 
and the statistical error bars on the points also apply to these curves. The redshift of Mrk 501 
z s = 0.034 was used in calculating the optical depth for each model. This figure is based on 
Figure 10 of [54]. 

is more absorption nearby in the Late SF model because the starlight in this model 
is less diluted by the expansion of the universe. It is possible that measuring the 
transparency of the universe to 7-rays at ~ 0.1 TeV with a number of sources at 
various redshifts can provide a strong probe of star formation, although there are 
uncertainties due to extinction by dust. 

V OUTLOOK 

The semi-analytic modelling of the EBL described here follows the evolution 
of galaxy formation in time. Forward modelling is a more physical approach 
than backward modelling (luminosity evolution). Pure luminosity evolution (e.g., 
[55-57]) assumes that the entire evolution of the luminosity of the universe arises 
from galaxies in the local universe just becoming brighter at higher redshift by some 




FIGURE 7. The 7-ray attenuation edge. The redshift where the optical depth reaches unity is 
shown as a function of 7-ray energy for each of the four ACDM models considered in Fig. 4. Also 
shown for the Kennicutt IMF is the redshift where the optical depth equals 0.5, 2, 3, and 5. 



power of (1 + 2) out to some maximum redshift. It effectively assumes that galaxies 
form at some high redshift and subsequently just evolve in luminosity in a simple 
way. This is at variance with hierarchical structure formation of the sort predicted 
by CDM-type models, which appears to be in better agreement with many sorts of 
observations. 

An alternative approach to modelling the EBL has been followed by Pei and 
collaborators [58-60], in which they find an overall fit to the global history of star 
formation subject to constraints from input data including the evolution of the 
amount of neutral hydrogen in damped Lya systems (DLAS). Their first attempt 
[58,59], which was used as the basis for EBL estimates by [35,61], was somewhat 
misled by the sharp drop in the DLAS hydrogen abundance from redshift z ~ 3 to 
z ~ 2 reported in [62]. With more complete data on DLAS (see, e.g., Fig. 14 of 
[63]) the z = 3 neutral hydrogen abundance is lower and almost constant from z = 2 
to 4. The latest paper by Pei et al. [60] takes a variety of recent data into account. 
Their approach is to follow the evolution of the total mass in stars, interstellar 
gas, and metals in a representative volume of the universe; they assume a Salpeter 
IMF. By contrast, the semi-analytic methods we use follow the evolution of many 
individual galaxies in the hierarchically merging halos of specific CDM models, 



here ACDM. Despite the differences in approach, and the fact that [60] assumed 
Q m = 1 and Hubble parameter h = 0.5, their results are broadly similar to those 
from the semi-analytic approach (see their §4.4; for our semi-analytic approach to 
modelling DLAS, see [64]). In particular, their EBL is similar to our old results [30] 
for the Salpeter IMF. Our EBL results presented here are higher in the near-IR and 
more consistent with the direct determinations [45,46]; they are also higher in the 
mid-IR, probably mainly because of the warm dust and PAH features in our dust 
emission model. It will be interesting to see whether further development of the 
global approach of Pei et al. and of the semi-analytic approach lead to convergent 
results. 

As our calculations show, the EBL, especially at £ 1 /im and ^ 10 /im, is sig- 
nificantly affected by the IMF and the absorption of starlight and its reradiation 
by dust, as well as by the underlying cosmology. The cosmological parameters 
are becoming increasingly well determined by other observations. As data become 
available on 7-ray emission and absorption from sources at various redshifts, es- 
pecially from the new generation of Atmospheric Cherenkov Telescopes and the 
new 7-ray satellites AGILE and GLAST, these data and their theoretical inter- 
pretation will help to answer fundamental questions concerning how and in what 
environments all the stars in the universe formed. 
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